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We consider tunneling of electromagnetic waves through a polariton band gap of a 1-D chain 
of atoms. We analytically demonstrate that a defect embedded in the structure gives rise to the 
resonance transmission at the frequency of a local polariton state associated with the defect. 



The optical properties of materials with band gaps in their electromagnetic spectrum have recently attracted, a 
great deal of attention. Iti-jwas suggested that fundamental electromagnetic processes such as spontaneous emissionJj'El 
photon-atom interaction,™] and optical energy transfero are strongly modified at band gap frequencies. Photonic 
crystals, which are periodic structures with a macroscopic period,ta present one of the primary examples of systems 
with electromagnetic band gaps. The periodicity in photonic crystals gives rise to allowed and forbidden bands for 
electromagnetic waves in basically the same manner as periodicity in the arrangement of atoms causes the band 
structure for electrons in solids. 

An important property of photonic crystals is an occurrence of local photon states with frequencies inside band 
gaps, when the periodic structure is locally distorted. The fact that an isolated defect in an otherwise perfect periodic 
crystal can give rise to local modes with frequencies in forbidden gaps of a host structure is well known in solid state 
physics. Local photons are similar in many aspects to other types of local states: Their frequencies always belong 
to forbidden gaps; in 3-D systems they split off the continuous spectrum only if the "strength" of a defect exceeds 
a certain threshold ;Bcj and by changing the type of the distortion one can control the position of the states inside 
the gap. It is essential, however, that while all other local states appear due to microscopic (of atomic dimensions) 
defects, local photons require both a macroscopic host structure and its macroscopic distortion. This fact is obviously 
due to the large wavelengths of electromagnetic waves in frequency regions of interest. 

Recently it was suggested that polariton gaps ia regular dielectrics with strong polar properties could generate effects 
similar to those attributed to photonic crystals.cliL3 For example, bounded photon-atom states originally considered 
in Ref. [ || for photonic crystals were proposed in Ref. [ |) for frequencies within a polariton gap. These states are 
formed by an optically active impurity atom, which possesses its own dipole moment different from that of the host 
atoms. j— - 

There exist a different type of local photon states, that is not associated with inner optical activity of impurities. E£l 
Such local states are analogous to defect modes in photonic crystals, and similar to them can be used to modify 
spectrum of radiation of the optically active impurities. At first glance it seems impossible, since electromagnetic 
waves would not interact with microscopic objects without inner optical activity. However, it was shown in Ref. [ 
|lO|j , that a regular microscopicdeiect embedded in a crystal lattice gives rise to local states with frequencies within 
the polariton gap, which are a mixture of electromagnetic component with excitations of a crystal responsible for 
the polariton gap. The most remarkable property of the local polaritons is absence of a threshold for localization in 
isotropic 3-D systems, which is due to a strong van Hove singularity in the polariton density of states at the gap edge, 
see details in Refs. [ [lO|JTl| . 

In this paper we show that the local states considered in Ref. [ |To| give rise to an interesting possibility of 
resonance tunneling of electromagnetic waves through a polariton bandgap. We would like—to emphasize that this 
tunneling proccess is remarkably different from both quantum mechanical electron tunnelingEij and photon tunneling 
through photonic crystalsJH The later proccesses can be characterized as the result of interaction between excitations 
and defects of comparable scales (electrons and impurities, electromagnetic waves and macroscopical distortion of 
photonic crystals). The tunneling studied in our paper occurs due to interaction between electromagnetic waves with 
macro scopically long wavelength and microscopical impurities. This process becomes possible owing to participation 
of the phonon component of polaritons, which mediate electromagnetic wave propagation. A similar role is played 
by excitons in resonance scattering of exciton-polaritons due to impurities with a short-range potential considered by 
Hopfield in Ref. [[[§. 

We present an exact analytical solution for the transmission coefficient of a scalar wave propagating through a l-D 
chain of noninteracting atoms containing a defect. These atoms are coupled to the wave due to a dipole moment 
caused by their mechanical vibrations. The spectrum of the coupled excitations of the chain and the field, polaritons, 
have a spectral gap where the excitations can exist only in an evanescent form. We show, however, that a defect, 
embedded in such a structure, results in the resonance tunneling of waves with the transmission coefficient independent 
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of the chain's length and being of the order of magnitude of one. One dimensional models usually describe tunneling 
processes fairly well because, by virtue of tunneling, the propagating wave is effectively confined in the transverse 
directions. In our particular situation it is also important that-tfae. local polariton states (transmitting centers) occur 
without a threshold in 3-D systems as well as in 1-D sy stems. Eiliil 

A similar model has been studied numerically in Ref. [ ^| , where a direct interaction between atoms of the chain 
(leading to the spatial dispersion of the chain's excitations) has been taken into account. The results of that paper 
suggest that though the spatial dispersion brings about some new features, it does not effect the existence of the 
resonance, justifying our neglect of the inter-atomic interaction. 

The atoms in our system are represented by their vibrational polarizability /3„, where subindex n represents the 
position of the atom in the chain. The polarizability is given by 
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where a is a coupling parameter between the dipoles and the field, and fl„ represents an atom's vibrational frequency. 
The defect in our model differs from host atoms in this parameter only, so = f2§ for all sites except of one occupied 
by the defect, where fl^ = fif. Polaritons arise as collective excitations of dipoles (polarization waves) coupled to the 
electromagnetic wave E(x n ) by means of a coupling parameter a. The electromagnetic subsystem is described by the 
following equation of motion 
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where the right hand side is a polarization density caused by atomic dipole moments, and c is the speed of the wave 
in vacuum. The coordinate x in eq. (|^) goes along the chain with the interatomic distance a. 

We first derive an equation for the frequency of the local polariton state in the 1-D situation. The one-dimensional 
nature of the model allows us to approach the problem microscopically and to take into account short-wave components 
of the field, including those beyond the first Brillouin band, exactly. Passing to the longwave limit at the last stage 
of the calculations we avoid non-physical divergencies and renormalization procedures of the kind used in Ref. [ Bj. 
The equation for an eigenfrequency of the local mode is 
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where Ail 2 = — £7q. It has a real- valued solution only if the frequency falls into the gap between the upper and 
the lower polariton branches. The integral in eq.@ can be calculated exactly to yield 

w = Of - <r , (4) 
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where we passed to the longwave limit toa/c <C 1, and introduced parameter d 2 — Ana/ a, which determines the width 
of the polariton gap between Oq and fig + d 2 . The second term in eq.(f|) is small for realistic values of the parameters, 
therefore, the frequency of the local mode is only slightly different from the defect frequency Of. As we shall see 
below, this fact has a deep impact upon the transmission frequency profile of the chain. 

The field in the eigenmode corresponding to the frequency determined by eq. (0) exponentially decreases away from 
the defect site: 

E = E def exp [-na(n - no)} , (5) 
where k is an inverse localization length of the state, which in the long wavelength approximation is given by 
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In order to consider transport properties of the model one has to subject eq.(Q) to the standard boundary conditions. 
We assume that incident and transmitted electromagnetic waves propagate in vacuum so that the transmission, t, 
and the reflection, r, coefficients are defined in the usual way 

dE dE 
E(0) = E in (l + r); — = ik u E in (l - r); E(L) = tE in exp (ik w L); — = ik w tE in exp (ik u L), (7) 
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where k u = lo/c is a wave number of the electromagnetic wave in vacuum, and L is the length of the chain. Our first 
goal in treating the problem of resonance tunneling is to convert the differential equation <M) into the discrete form. 
We can do so considering separately free propagation of electromagnetic waves between sites and scattering due to 
the interaction with a dipole moment at the site. Let E n and E' n be the magnitude of the electromagnetic field and its 
derivative just after scattering at the nth site. The electric field E remains continuous at a scattering site, while its 
derivative undergoes a jump, which is equal to — 47r/c 2 P n . Finally, one can derive the system of difference equations 
that can be written with the use of the transfer matrix, T, in the form: 

v n+ i = T n v n , (8) 



where we introduce the column vector, v ni with components E n , D n (D n — E^/k^). The transfer matrix, T, 
describes the propagation of the vector between adjacent sites: 
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The dynamical state of the system at the right end of the chain, which is represented by the vector vn, can be found 
from the initial state at the left end, vq, by means of the repetitive use of the transfer matrix: vjy — Y\± T n v$. In the 
case of a single impurity, the matrix product is conveniently calculated in the basis, where the matrix T for a host site 
is diagonal. With the use of the boundary conditions given by eq. (0) one finally arrives at the following expression 
for the transmission coefficient, i, in the limit kL > 1: 

t = *° (10) 

(1 + e) + i exp (-ik u L) Tt cosh [kcl(N - 2n + 1] ' 



where to is the transmission coefficient in the pure system, T = 47rfc w £/3/(sin(fc w a)v A D), D = (47rfc w /3) 2 + 
167rfc w ,3cot(a/c w ) — 4, and e — Ank^ (fidef — 0) /"vD. The last parameter reflects the difference between host atoms 
and the impurity, and is equal to 
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The transmission coefficient to of the pure system is given by 

t = 2e^exp(-ttL) 

l-^(2-47rfc w /3cotfc w a) 

and exhibites a regular exponential decay. 

Eq.([l0|) describes the resonance tunneling of the electromagnetic waves through the chain with the defect. The 
resonance occurs when 

l + e = 0, (13) 

with the transmission becoming independent of the system size. Substituting the definition of the parameter e from 
eq.([ll]) in eq.(|l3|) one arrives at the equation identical to eq.(|J) for the frequency of the local polariton state. Typically, 
for resonant tunneling, the transmission takes a maximum value when the defect is located in the middle of the chain, 
N — 2«o + 1 = 0, and in this case 
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The width of the resonance is proportional to Tto and exponentially decreases with an increase of the system's size. 
In the longwave limit ak u <C 1, eq.(|l4|) can be rewritten in the following form 

2 ' '^-2^-) , (15) 
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where ui r is the resonance frequency satisfying eq. ( |l3| ) . It is interesting to note that the transmission becomes exactly 
equal to one if the resonance frequency corresponds exactly to the center of the polariton gap. This fact has a simple 
physical explanation. For cj 2 = + <i 2 /2 the inverse localization length k becomes exactly equal to the wave number 
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Lo r / c of the incoming radiation. Owing to this fact the field and its derivative inside the chain exactly match the field 
and the derivative of the incoming field as though the optical properties of the chain are identical to those in vacuum. 
Consequently, the field propagates through the chain without reflection. 

Equations (|l(j) and ( |l4|) demonstrate that the frequency profile of the resonance is considerably different from a 
corresponding Lorentztian profile for electronic and other known instances of tunneling. This occurs because the 
parameter e diverges at to — £l\ causing the transmission to vanish. At the same time the resonance frequency oj r is 
very close to fii as it follows from eq.(H). This results in a frequency dependence of the transmission strongly skewed 
toward lower frequencies. 

Having solved the transmission problem, we can find the magnitude of the field inside the chain in terms of the 
incident amplitude, Ei n , at the resonance frequency. Matching the field in the local polariton state given by eq. (|^) 
with the outcoming field we have for the field amplitude at the defect atom 

Ed = Ei n t exp(— i/cL) exp [(N — n^na] . (16) 

For \t\ being of the order of one in the resonance, this expression describes the drastic exponential enhancement of the 
incident amplitude at the defect side due to the effect of the resopance tunneling. This effect is an electromagnetic 
analogue of charge accumulation in the case of electron tunneling.113 

Resonance tunneling is very sensitive to the presence of relaxation, which phcnomcnologically can be accounted for 
by adding 2i^to to the denominator of the polarizability f3. This will make the parameter e complex valued, and the 
resonance condition Re(e) = — 1 may only be fulfilled if the relaxation is as small as 7 < (ad 2 )/ (Ac). This inequality 
has a simple physical meaning: it ensures that the distance between the resonance frequency and f2i, where the trans- 
mission goes to zero, is greater than the relaxation parameter 7. This is a strict condition that can only be satisfied for 
high frequency oscillations with large oscillator strength in molecular crystals with large molecules in an elementary 
cell, and respectively large values of the interatomic spacing a. Another interesting opportunity can arise in so called 
atomic optical lattices, where atoms, trapped by a laser beam, form a lattice with spacing practically equal to the 
wavelength of the trapping field.H3 However, taking into account the spatial dispersion can lead to a more favorable 
situation for the tunneling resonance in our model. Numerical results of Ref. [ |l2j show that spatial dispersion does 
not change the transmission properties significantly. Therefore, one can rely upon eq.([To|) to estimate the effect of the 
dissipation in the presence of the spatial dispersion, assuming that it only modifies the parameter e . Eigenfrequency 
of the local mode in the presence of the direct inter-atomic interaction is calculated exactly. The interaction moves the 
resonance frequency farther away from fii undermining the influence of the damping. Adjusting e correspondingly one 
can arrive at the following new condition for the resonance to survive the relaxation: (7^1)/$ < 1, where $ is the pa- 
rameter of the inter-atomic interaction; it can be estimated as a bandwith of the polarization waves (in terms of squared 
frequencies). This condition can be fulfilled, even for phonons with a relatively small negative spatial dispersion, and 
becomes even less restrictive in the case of Frenkel excitons in molecular crystals. The imaginary part of e will prevent 
the exponential factor to in eq.(|lfj|) from canceling out at the resonance. This restricts the length of the system in 
which the resonance can occur. The requirement that Im(e) be much smaller than leads to the following condition 
L <^1/(k\ ln[Im(e)] |), with Im(e) being again of the order of min[(47c)/(ad 2 ), (7^1)/$]. 

Concluding, we showed that a regular defect atom without internal degrees of freedom and an optical activity results 
in resonance tunneling of electromagnetic waves through a polariton gap . Though we considered the one-dimensional 
model one can expect that the existence of the effect does not depend upon dimensionality because tunneling transport 
is virtually one-dimensional, and the polariton local states are present in the system of any dimension. 
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